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We discuss a non-equilibrium statistical system on a graph or network. Identical particles are 
injected, interact with each other, traverse, and leave the graph in a stochastic manner described in 
terms of Poisson rates, possibly dependent on time and instantaneous occupation numbers at the 
nodes of the graph. We show that under the assumption of constancy of the relative rates, the system 
demonstrates a profound statistical symmetry, resulting in geometric universality of the statistics 
of the particle currents. This phenomenon applies broadly to many man-made and natural open 
stochastic systems, such as queuing of packages over the internet, transport of electrons and quasi- 
particles in mesoscopic systems, and chains of reactions in bio-chemical networks . We illustrate the 
utility of our general approach using two enabling examples from the two latter disciplines. 

Analyzing the statistics of currents generated in an open driven system that consists of multiple degrees of freedom is 
the holy grail of non-equilibrium statistical physics. The problem typically emerges in the context of multiple particles 
traversing a medium, subject to stochastic noise, and also interacting with each other. To study an elementary 
current means to count the number of particles passing through a respective line element of the medium during 
a fixed observation time. Repeated a multiple number of times, these measurements of the current show variable 
results due to the intrinsic stochastic nature of the underlying processes. To quantify these statistical variations one 
studies the Probability Distribution Function (PDF) of the current. It may also be useful and instructive to study 
simultaneously multiple currents, even the complete set of currents, related to all linear elements of the medium. When 
the observation time is sufficiently long, in particular much longer than a typical time for a particle to traverse the 
system, one expects the PDF of the currents to obey a self-similar, so-called Large Deviation (LD) form, expressing the 
fact that, the longer the observation time is, the smaller the fluctuations in current density (total current accumulated 
over the observation time divided by the latter) are. The LD is the most general, yet not the only universal statement 
one can make about the PDF of currents. Thus, one popular formulation of the so-called Fluctuation Theorem (FT) 
suggests that under rather general conditions the PDF of directional currents (which can be positive or negative, as 
counting the particles traversing a line element in different directions with different signs) also shows an additional 
profound symmetry, e.g. relating to each other the statistics of large negative and large positive currents [TrE] . In this 
manuscript, we suggest a set of even more general (than FT) statements about the statistics of currents, observed in 
a variety of open many-body (i.e. consisting of multiple interacting particles) stochastic systems. 

Examples of the systems, in which our consideration applies, are wide-spread in physics, biology, and operations 
research. Fig. [I] shows a graph of transitions between the nodes traversed by multiple jobs, cars, internet packets, or 
particles, which are processed/delayed at the nodes, e.g. competing for the servers' time, with the delay representing 
interactions for this example from Queuing Theory/ Operations Research [5H3]. In statistical physics similar models 
are called zero-range processes [TOl HT] . Bio-chemical networks represent another area of application for the current 
(counting) statistics [T2H25] . The multiplicity of reactions for a given bio-chemical setting can be represented in 
terms of a graph of transitions between different bio-chemical states, e.g. as illustrated in Figs. [2j Even when 
the graph is deduced reliably, the transition (kinetic) rates often remain hidden, subject to inference/reconstruction 
from macroscopic measurements of counting statistics, which have recently extended the arsenal of experimental 
techniques available in the field. Yet another application domain which requires very precise evaluation of statistics 
of currents belongs to the field of mesoscopic electronics, e.g. electric circuits in systems made of quantum dots and 
tunnel junctions [2rJH251 [251 1251 [25ll33| . with an example illustrated in Fig. [3j Better understanding of the electron 
transport/current in such devices is instrumental for designing fast and reliable information processing systems for the 
future. These applications and examples will help us to illustrate the power of our main result, stated and explained 
in the next Section. 



I. MAIN RESULT: GEOMETRIC UNIVERSALITY OF CURRENTS 

In this manuscript we consider identical particles traversing the system, described by a connected open directed 
graph, S = (So, Si), for example shown in Fig. [I] where So and Si denote the sets of the graph vertices and edges, 
respectively. Particles enter the system from the exterior reservoir of infinite capacity, where Aoi is a time-independent 
Poisson particle injection rate at node i. Upon arrival at an entry node the particle advances by jumping between the 
nodes along the directed edges, each characterized by the transition rate, = 7ijAifc(f; n(t)), with Xik being the rate 
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FIG. 1: Example of an open network of interacting particles represented by a directed graph. The sample graph consists of 
four vortices/stations, labeled 1,2,3,4, with label reserved for an external (reservoir) node. Transitions between the nodes 
are shown as directed edges. Loops (self-loops, as 1 — > 1) are allowed. Each graph edge (i, k) with i 7^ is equipped with 
a transition rate Xik, describing the rate of particle departure from site i per current occupation number rii on the site. In 
this Letter we consider the most general case of Xik{n(t);t) being finite, and otherwise depending arbitrarily on time and the 
instantaneous global state n(t) of the system. Aoi stands for the time-independent rate of injection from the reservoir. We also 
allow two or more directed edges to connect the same vertex. In this context we mark the edges by distinct subindexes, like 
edges 23 1 and 23 2 in the Figure. 

per particle that occupies node i. The internode per-particle rate, X^, is finite and may depend on time explicitly, as 
well as implicitly via the dependence on the system state n = (rii\i £ 5o), with rii(t) being the number of particles 
at node i at time t. We assume no intrinsic limitation on the number of particles residing at any node in the system, 
i.e. rii may be any nonnegative integer. A particle may leave the system along any of the outgoing links, with the 
corresponding departure rates Aj = TiiXio, the per-node rate A^o being also generally dependent on both n{t) and t. 
Such a stochastic system of many identical interacting particles is described via the following Master Equation (ME) 
for the instantaneous distribution function of n: 



where n + i ~ (■ ■ • , + 1, - • • ), n_i = (• • ■ , ni — 1, • • • ) and n + i ; _fe = (• • ■ ,m + 1, • • • , n& — 1, • • • ). We assume that the 
graph and the rates are chosen so that the solution of Eq. ([!]) converges to a non-singular quasi-steady distribution of n 
in some finite time, r, which can be also viewed as the typical time a particle spends in the system. To characterize the 
non-equilibrium features of this steady state one may also want to know the statistics of the vector J — ( Jik\(i, k) £ Si) 
of currents, with the edge-labeled components, Jik, defined as the number of particles, that traverse through the edge 
(i, k) during time T. For the observation time, T, that is long compared to the correlation time, t, the PDF of the 
currents obeys the following Large Deviation (LD) form: - ln(T(J|T))/T = S(J /T) + o(t/T), where §(J/T), often 
referred to as the LD function, is a convex function of its multi- variant argument (in some physics and mathematical 
physics literature it may be also referred to as the Cramer function). 

The main statement of this manuscript, which we call the Geometric Universality of Currents (GUC), suggests 
that for T»t and under the condition that the tested currents are not very large (so that the average n conditioned 
to the value of currents is non-singular and finite), the LD function is invariant with respect to any (e.g. time and 
n(t) dependent) transformations that keep the graph intact and the branching ratios 
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constant, i.e. the transition probabilities 9ik are independent of the current time t and the current state n(t) of the 
network. In particular, the statement means that §( J /t) can be calculated by considering an auxiliary problem with 
Xik(t; n) — > Oik, thus replacing the interacting system with its much simpler non-interacting surrogate. This reduction 
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makes the evaluation of the LD function computationally easy. We label the universality, geometric, to emphasize 
its interpretation in terms of the statistics of open single-particle trajectories, and specifically of their geometries. 
Some geometric aspects of this analysis were already discussed in [31] for a particular case of a queuing model with 
time-independent rates. 

Our proposed GUC also allows generalization to the case in which the injection rates dependend on time arbitrarily, 
for example in a quasi-periodic fashion, so that the time averaged rates Xoi = T" 1 dt\oi(t) are finite. Then, under 
the same conditions [universality and independence of time for all #'s in Eq. ([2]), and finiteness of the queues] as for 
the above statement of GUC, we find that statistics of currents is universal with the LD function calculated, as if 
the injection rates were constant and equal to Ao^. We refer to this extension of GUC as the GUC for dynamic 
pumping. 

In what follows we will: (a) Discuss a set of examples from the literature, e.g. on queuing networks, networks of 
bio-chemical reactions, and electron transport over mesoscopic networks, in which many-body setting and associated 
statistics of currents are of interest, (b) Sketch a proof of the GUC. (c) Illustrate applications of GUC to bio- 
chemical networks and systems of quantum dots, (d) Conclude the manuscript by describing possible extensions and 
generalizations of GUC to other characteristics (such as the joint distribution functions of occupations and currents) 
and more general settings. 

II. MODELS 

We begin this Section by explaining a basic setting in the Operation Research (Queuing Networks) and then proceed 
with a discussion of our two enabling applications: one from the field of bio-chemical networks, and the other being 
related to electron currents in mesoscopic systems of quantum dots. 

Queuing networks. Introduced and studied in the work of Jackson [6] and Kelly [7], this subject became classical 
in Operations Research (OR), see e.g. the textbooks [E1IH]- Even though stated in somewhat different terms, the 
basic equation of the field, concerned with maintaining finite multiple queues in industrial conveyors, calling centers, 
transportation, and internet systems, is still the ME ([I]) with a specific choice of the rate dependence on the vector n 
of queues. In particular, the so-called Jackson network model assumes \k to be equal to a time-independent constant 
for ni < m, and zero otherwise, with m, interpreted as the number of servers available at node i. The particles (jobs, 
calls, cars, or packets) compete for the servers that may be busy serving other particles and thus interact, unless 
the uii are infinite at all nodes of the network. In the latter case the particles are "free", i.e. they do not interact. 
The queuing network is called stable if a statistically steady solution of ME ([lj exists, with all queues being finite. 
It was realized quite early in the history of this field [7 that if a network is stable, the stationary PDF of queues 
becomes factorized into the so-called product form, P s ^(n) = Y\i Pi{ n i), where pi are the marginal distributions at the 
nodes. A similar decomposition has been also discovered by Spitzer |10j in the related mathematical physics context 
of interacting Markov processes. Statistics of currents in Jackson networks was studied much less in the classical OR 
literature. See e.g. [5J[35]- Thus, the LD function of currents in the (static) Jackson networks was analyzed only 
recently |34j . In this regards, this manuscript extends the approach and technique developed in |34| to a more general 
setting, in which transition rates in the ME ([I]) depend not only on the occupation number of the outgoing vertex, 
but also on the global state of the system, characterized by n(t), and may also depend explicitly on time. 

Biochemical Networks. When the number of molecules involved in a bio-chemical reaction is not too large, e.g. 
when it is of the order or less then 100, the effects of noise on reaction kinetics of a moderate-size system (often called 
a mesoscopic system) are significant. For example in vivo reactions inside living bacteria [16 18, 36 with hidden 
kinetics (reaction rates are not known) is an exemplary mesoscopic system of this type. Effects of fluctuations in 
such systems are often analyzed in vitro using single-molecule experiments and methods of fluorescence correlation 
spectroscopy [19 25 . These and other related techniques try to reconstruct the complex stochastic mechanisms 
of biochemical networks. In particular, the techniques are capable of uncovering new microscopic information on 
reaction kinetics, inaccessible via standard macroscopic (bulk) experi ments. This is achieved by measuring variances 
and autocorrelation functions of the reaction events. Examples of observations, which became available with the 
invention of single-molecule techniques, include the discovery of new (not known before) internal states of enzymes 
[22j . The opportunities provided by these techniques are exciting, however they are also accompanied by significant 
challenges in interpreting the statistical measurements. Direct numerical simulations of even simple biochemical 
processes are prohibitively expensive and inconclusive for reconstructing full statistics of reaction events [37], thus 
emphasizing the importance of alternative analytical methods. This manuscript contributes to the task of developing 
analytical methods to describe delicate and important fluctuation effects in mesoscopic bio-chemical reactions. We 
specifically consider the following exemplary model. 

Model 1: Noise in Enzymatic Reactions: Consider the conversion of one type of molecules, called the substrate, 
into another type, called the product, via a complex Michaelis-Menten (MM) enzymatic reaction [38] ■ This process 
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FIG. 2: (a) Set of reactions that include a complex Michaelis-Menten process, as well as enzyme creation and degradation, (b) 
Equivalent Jackson network with kinetic rates which are proportional to the number of molecules, represented by the nodes. 
First and second nodes represent free enzymes and enzyme-substrate complexes, respectively. 




FIG. 3: Electric circuit made of three quantum dots (QD), controlled by gate voltages, Vi, (i = 1,2,3). Quantum dots are 
coupled to each other and to the leads via tunnel junctions with barriers, characterized by the parameters, gij = g^ and gt, 
respectively. Leads (green) are characterized by their chemical potentials, Iinset: Experimental realization of a double 
quantum dot from Fig. 7(a) in Ref. [31] • Electron transport between source (S) and drain (D) proceeds through tunneling via 
two dots with potentials, tunable by gate voltages, Gi and G r - The barrier strengths at the tunnel junctions between two dots 
and between the dots and the leads S and D are controlled by additional gate voltages Gl, Gr and Gc- 



involves interactions of a substrate molecule with another molecule, called the enzyme, that results in the creation of 
an enzyme-substrate complex. Let S, P, E, and C denote the substrate, the product, the enzyme, and the enzyme- 
substrate complex molecules, respectively. The complex, C, is created from E and S, followed by splitting into either 
E + S or E + P. Assume the latter reaction is observable. This is possible, e.g. via attaching the so-called Green 
Fluorescent Protein (GFP) tag that fluoresces each time a product molecule is created [22J. The light intensity is 
proportional to the number of product molecules. Measurements of the average intensity, its variance, and higher 
cumulants provides valuable information on the statistics of the number of transitions in the C — > E+P sub-chain. The 
enzyme molecules are also permanently generated and undergo degradation. To summarize, the full set of elementary 
reactions that characterize the model is shown in Fig. [2ja). It represents, (1) creation of enzyme molecules, — > E, 
with rate fc + ; (2) degradation of the enzyme-substrate complex C — > with rate fcc[C]; (3) degradation of enzyme 
molecules, E — > 0, with rate };e[E]; (4) Conversion of S into P via a Michaelis Menten reaction, which involves three 
sub-processes: (a) Creation of an enzyme-substrate complex E + S — >• C with rate fcif.EjfS 1 ]; (b) Reverse reaction of a 
complex, splitting into free enzyme and substrate molecules, C — » E + S 1 , with rate A:_i[C]; (c) Irreversible splitting 
of a complex into free enzyme and a product molecules, C — > E + P, with rate fc 2 [C]. Here [. . . ] represents a standard 
bio-chemical notation for abundance [37], i.e. the number of molecules of the given type in the system. We assume 
that the substrate molecules are supplied in macroscopic quantities, thus keeping the corresponding concentration 
[S] = const. Note that the described stochastic model of bio-chemical kinetics can be viewed in terms of stochastic 
transitions in the reaction graph shown in Fig. [2^b) , where the number of transitions through a given link represents 
the number of reaction events during the observation time. Assuming all transitions to be Poisson, we observe that 
the set of bio-chemical reactions in Fig. 2ja) can be viewed as an instance of a general stochastic network model 
with the property (J2|. Specifically, the relations between the parameters of the bio-chemical models with fluctuating 
numbers of enzymes and the notation for the rates in our general network model (see e.g. Eq. ([I])) are as follows: 
A i = k+, Aio = k E [E], A 02 = 0, A 20 = k c [C\, Ai 2 = h[E][S], A 21 i = fc_i[C], A 21 2 = k 2 [C]. 

Transport in Driven Mesoscopic Systems with Detailed Balance. Many mesoscopic systems in physics and 
chemistry can be modeled in terms of stochastic transitions that satisfy the Detailed Balance (DB) constraints with 
the participating particles entering and leaving the system via a reservoir, or a system of reservoirs. DB guarantees 



FIG. 4: "Queuing-style" (oriented) graph representation for a system of three quantum dots shown in Fig. P 



that for time-independent kinetic rates the system attains the equilibrium Boltzmann-Gibbs distribution. In a system 
with DB, an allowed transition from state i to state j is always accompanied by the reversed transition. A transition 
between two states i and j can be viewed in terms of a barrier, with the transition strength characterized by a 
symmetric characteristic gij = gji, independent of the system state n(t) (number of particles at a given node at a 
given time). Nodes are represented by the local trap potentials with energies, Ei and Ej, that control the asymmetry 
between the forward and backward transitions. Energies may generally depend on the state (which represents multi- 
particle interactions) and time. It is convenient to adopt the Arrhenius parameterization of the DB per-particle kinetic 
rates transitions Xij — Sigij and Xji = e^g^ for the and transitions, respectively, with Ei — e@ Ei , and /3 
being the inverse temperature in energy units. We also assume that all nodes of the open system are exchanging 
particles with the (external) reservoir. In fact, it may be convenient to introduce individual reservoirs for the nodes of 
the "state" graph, and characterize them in terms of their chemical potentials, fj,i, and the barrier strengths, gi. By 
analogy with the inter-node transitions we assume gt to be time- and state- independent, while the chemical potentials, 
fii, may depend on time. Kinetic rates for transitions from node i to its reservoir and vice versa are A^o = £;<7i and 
Aoi = e i<?i! respectively, with = e^ 4 '. 

For systems with DB, described above, the symmetry property ([2| is naturally built in, provided the barrier 
parameters g^ and gi are kept time- and state-independent. Stated differently, the barrier parameters are single- 
particle, i.e. non-interacting, characteristics. Many-body interactions are represented by a (possibly complex and 
generally arbitrary) dependence of the node activation rates £ = (ei, S2, ■ ■ ■ , £jv)i on the system state n(t). Moreover, 
the non-equilibrium/driven nature of the system state originates from additional explicit time dependence of the 
activation rates e = (ei, £2, . . . , eat) and e = (ex, 62, ■ . ■ ejy) for the nodes and reservoirs, respectively. Periodic or 
quasi-periodic driving are two possible and representative examples. Measurements of the full counting statistics of 
currents in driven nanoscale devices, with otherwise DB conditions on kinetic rates, represent a feasible state-of-the-art 
in mesoscopic experiments in electronics [321 1331 140] and soft matter |41j . 

Model 2: Transport in a network with DB is implemented in an electric circuit, such as in Fig. [3} that includes a 
number of semiconductor quantum dots populated by electrons which are allowed to hop through the tunnel junctions 
between the quantum dots and corresponding reservoirs. A multi-reservoir queuing-style graph representation that 
corresponds to the electric circuit in Fig. [3] is shown in Fig. [4j Electric circuits made of several quantum dots with 
tunable parameters have been implemented experimentally [28-30J. Such systems are used to create new devices 
capable of few-electrons information processing [55], as well as showing useful features of quantum dot ratchets [29] . 

Our main focus for Model 2 is the possibility of observing the so-called pump effect that refers to efficient generation 
of currents induced by external periodic changes of the parameters, such as gate voltages that control the charging 
energies, Ei, in the quantum dots [42]. The pump effect can be used to build a device, called an electronic turnstile 
[30j . that involves a circuit with quantum dots, driven by periodic gate voltages at zero bias between the leads. When 
the gate voltages are fixed, the electric current is not generated, however, under proper conditions, periodic driving 
generates currents, whose magnitudes can be controlled with high precision, thanks to the experimental ability to 
control the number of electrons pumped per cycle of the driving protocol |43) . We consider such periodically driven 
electronic circuits in the regime of sufficiently high temperatures, so that the number of available states inside a 
quantum dot is large compared to a typical number of electrons in a dot. In this regime quantum coherence effects 
and exclusion interactions due to Pauli principle can be disregarded, so that electronic transitions through tunnel 
junctions are thermally activated and the corresponding rates can be described using the Arrhenius parametrization. 



FIG. 5: Two geometrically distinct paths for single-particle trajectories. 



Kinetic rates are controlled by application of gate voltages to quantum dots, as shown in Fig. [3] Note, that our 
forthcoming discussion of currents in the model accounts for many-body effects, originating, e.g., from Coulomb 
electrostatic interactions. 



III. MATERIALS AND METHODS 

This Section has two parts. We first give a sketch of the GUC proof and then illustrate the utility of our general result 
for two enabling applications, related to bio-chemical networks, and electron transport in quantum dots, respectively. 

A. Proof of GUC 

Our main result, GUC, can be rationalized in two different ways. The Lagrangian (Dynamic) approach is discussed 
below, while the complementary Eulerian (Static) consideration and the proof of GUC generalization to Dynamic 
Pumping are given in Section I of SI. 

1. Lagrangian (Dynamic) Approach 

Although in the interacting case the particle trajectories are strongly correlated, the main focus of our Lagrangian 
consideration will be on a single-particle trajectory, x = (s,t), fully characterized by the following finite sequences, 
s = (si,...,Sfe) and t = (to, t\ . . . ,tfc), where si and i;, with 1 < I < k, stand for position (node s/ G So) and 
time stamp of the particle leaving the station s; for the next station s; +1 along the path. We assume that to and 
tk are the times when the particle enters and leaves the system, respectively; the temporal stamps are time-ordered: 
< to < ■ ■ ■ < tk < T, where T ^ r is a sufficiently long time and r is the correlation time (estimated as the typical 
time a particle spends inside the system); a particle can visit the same station a multiple number of times; and the 
total number, k, of jumps can be an arbitrary integer. It will be important for the forthcoming discussion to separate 
the geometrical, s — (s 1; • • • , s^), and temporal, t — (i , • ■ • , ifc) characteristics of the trajectory. In particular, the 
total current generated by a single-particle trajectory, whose individual (per edge) component measures the number 
of times the trajectory passes through any given edge, depends on its geometrical component, s, only, and the current 
generated by a stochastic realization of our system is the sum of the single-trajectory currents. Thus, to describe a 
single-particle trajectory, for example one of the two shown in Fig. [5] we naturally pose the following question: what 
is the probability for a particle observed entering and leaving the network over sufficiently long time horizon [0; T] 
to choose this particular path, s? This probability should depend on the factors 9 SlSl+1 only and, therefore, it shows 
no correlations with other trajectories. The information contained in this probability distribution, T(s) is geometric 
and independent of the temporal details t of the trajectory. The assumed invariance (§ of our general model, with 
respect to any changes in the system keeping the transition (branching) probabilities On constant, guarantees that 
the probability of a path, CP(s) = Ilf=i^s!s i+ i (with s; +1 = 0), is universal/invariant as well. The probability of a 
temporal realization of a many-particle trajectory shows strong correlations between single-particle trajectories and 
can be explicitly expressed in terms of k Si (t; n(t)) = X^es^ 31 however, the generated currents are not sensitive 
to this complex temporal structure. 

We now consider the probability of observing a given vector of the current for the observation time, T (the component 
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of the vector associated with an edge counts the number of particles that have crossed the edge in time T.) 
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where ^(sj, • ■ ■ , Sg) is the probability of n particles being injected into the network at the corresponding nodes and 
Xife(s) simply counts the number of times the particles traversing the path s went through the edge (i, k). Obviously, 
Eq. ([3| is evaluated as an expectation over all possible single-particle paths. [Eq. ([3| is an asymptotic expression, 
valid for T ^ r only, because it completely ignores the trajectories, trapped in the system for periods longer than T]. 
Therefore, one concludes that asymptotically y{J\T) is a single-particle object which is also universal/invariant with 
respect to Eq. d2| in view of its linear relation to T(s) and the previously discussed universality of the latter. | 



Here we discuss two applications of GUC to Model 1 and Model 2. In the context of these two models GUC can be 
viewed as a no-go restriction on the ability to interpret experimental data on the statistics of reaction events (Model 
1), or to control electrical currents in networks with DB (Model 2). These restrictions follow from the fact that 
current statistics actually depends only on the time-independent relative probabilities, 6ij defined in Eq. |2| as ratios 
of transition rates, rather than the complete transition rates themselves. 

Model 1. Suppose that one performs measurements of the number of product molecules converted from the 
substrate by enzyme molecules via a set of enzymatic reactions. Assume that this experiment can be performed in 
a single-molecule setting and repeated to determine not only the average rate of the product creation, but also to 
measure and quantify mesoscopic fluctuations of the currents. One can attempt to use the measured data to extract 
information on the reaction mechanisms, e.g. to determine the kinetic rates for elementary reactions, involved in the 
enzymatic reaction. 

A version of the general GUC relevant for this model is: For steady (time-independent) rates of enzymes' creation 
and degradation, the counting statistics of the enzymatic reaction events does not provide full quantitative information 
on all kinetic rates of elementary reactions involved in the process. The number of independent parameters that cannot 
be determined is, at least, equal to the number of different states of the enzyme and the enzyme-substrate complex. 
This restriction follows from the observation that if we multiply all outgoing rates from any node of the graph of 
transformations by a constant factor, the full counting statistics of transitions through any link in the network will 
not change. Therefore, for each network node there is a set of parameters that cannot be extracted from the counting 
statistics. This general result is illustrated in Section B of SI using an example of the open-system Michaelis-Menten 
(MM) mod el of Fig. [2] The main no-go (prohibitive) conclusion of the open-system MM analysis is that the MM 
curve for the average current (along the upper C — > E link from Fig. [2| , presented as a function of [S] , depends only 
on two independent quantities, namely k 2 k + /k_ and the MM constant Kmm — (&2 + Measurements of the 

Fano factor (defined as the ratio of the second-to-first cummulants of the current) provide additional information on 
the ratio, k 2 /k^, thus allowing the value of k + to be extracted from a combination of the two measurements. However, 
the other two rates of the complex, k\ and cannot be determined even if fc_ is known independently from some 
complimentary experiment. 

Model 2: For a system with DB, it is convenient to consider undirected links that can be traversed in both 
directions, so that the current associated with an original (undirected) link is the difference between particles moving 
in the main (positive) and reverse (negative) directions. These currents will be referred to as bi-directional currents. 
Obviously, statistics of bi-directional currents are obtained from A(q) by setting q^ — 1/qji in the argument of the 
generating function. Then, according to GUC for dynamic pumping, the statistics for the vector of bi-directional 
currents, built of the components associated with undirected edges (we will still use q for the argument of the 
generating function, not to complicate the notation), is characterized by 



B. Applications of GUC 
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Eq. Q shows that in spite of the many-body interaction and explicit time dependence of the parameters e and e, the 
statistics of currents appears to be the same as in the auxiliary time-independent and non-interacting system with e 
expressing the effective constant chemical potentials of the reservoirs. We arrive at the following general (as applied to 
any open DB system driven externally in a periodic or quasi-periodic fashion) No-pumping Fluctuation Theorem 
(NPFT) relating currents in the original and reduced systems: If the barrier characteristics gij of the transitions, and 
gi of the reservoirs remain time-independent, while vectors e and e of the activation rates are driven arbitrarily in time 
(periodically or quasi-periodically) and e generally depends on the local node populations, then, provided the observation 
time is long compared to all correlation times in the system, the counting statistics of currents is indistinguishable 
from its counterpart in a reference non-interacting, £ = 1, system that stays in its steady state (not driven) with 
e being for the vector of the reservoir activation rates. This general statement contains no-pumping in its name, 
since it prohibits generation of currents that are different from the steady state currents under the aforementioned 
conditions for periodic pumping. In particular, NPFT guarantees that if the vector e of the reservoir activation rates 
has all the same entries, then all odd current cumulants remain zero, i.e. the current statistics are indistinguishable 
from their counterparts for the currents in thermodynamic equilibrium, even if the driving protocol explicitly breaks 
time-reversal symmetry. 

Extended discussions of the relation between NPFT and the previously known FT, as well as further illustrations of 
how NPFT can be used to build an efficient electronic turnstile ,30], are given in Sections II and III of SI, respectively. 

IV. CONCLUSIONS AND PATH FORWARD 

In this manuscript we have formulated and demonstrated a strong symmetry, referred to as Geometric Universality 
of Currents. GUC imposes restrictions on stochastic of currents generated over long times in the systems of interacting 
particles in open networks. We have shown that the GUC rests on the geometrical nature of currents, and is limited to 
large but not extreme values of the currents, so that almost all particles spend only a finite time within the network. 
We have observed that GUC extends only to the LD function (the leading, exponential term in the PDF) while the 
pre-exponcntial factor is not universal. 

GUC imposes strong restrictions on the amount of information one can extract by measuring the statistics of 
currents. From this viewpoint, studying the conditions under which GUC can be violated and extending GUC to 
more general settings can be expected to have increasingly important role in the future. This may take place via 
finding the dead ends and opening new venues for interpretations of experimental measurements in terms of the 
underlying stochastic mechanisms and microscopic parameters, as well as identifying efficient strategies for design 
and optimization of nanoscale electronic and biochemical devices. Therefore, we would like to conclude with a brief 
discussion of possible future generalizations. 

The most general setting, considered in this manuscript still has three basic restrictions: the observed currents 
are not too large, the number of particles occupying a node is unrestricted, and particles are identical, of only one 
kind. It has been recently shown [34] in the context of Jackson networks that even when the network is stable (i.e. 
when particles do not accumulate in the steady state), generation of atypically large currents is accompanied by 
accumulation of particles in the system. However, even in this regime of extreme currents, the LD form of the current 
distribution holds, CP(J|T) ~ exp(— TS( J/T)), and the LD function §(j) remains well-defined. The most probable 
way to reach these atypically large values of j is associated with accumulation of particles at some, so-called saturated, 
nodes. In this regime of extreme currents a finite fraction of particles will not leave the system during the observa 
tion time (in spite of the fact that the latter is much larger than the correlation time r), thus violating one of the 
key assumptions used to derive GUC. However, and as shown in (34] . the saturated nodes start acting as additional 
Poisson sources (reservoirs) , thus leading to survival of GUC in a partial form - with respect to variations in rates at 
the non-saturated nodes. To study this effect of the partial breakdown/survival of GUC under extreme currents is an 
interesting task for the future. 

Introducing limitations on the number of particles that can occupy the same node, as well as analyzing several types 
of particles that interact with each other differently at different nodes (the two regimes are known in the queuing 
literature as one of the finite waiting room, and one characterized by graph-inhomogeneous priorities, respectively), 
makes the analysis of currents much more complicated. Addressing these issues in the context of GUC should become 
yet another challenge for the future. Finally, let us mention that the GUC analysis may also be extended evaluating 
LD functions of currents in various interesting non-equilibrium cases quantitatively, e.g. in the spirit of (341 1441 [45] . 
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This supplementary material presents further details on the various results discussed in the main body of the paper, 
following the plan described hereafter. 

Section A] achieves two goals. First it describes in Section |A l| a complementary Eulerian (static) proof of the main 
result of the paper, Geometric Universality of Currents (GUC), and second it extends the Eulerian approach to prove 
GUC for dynamic pumping, which is done in Section [A 2| Section [B] describes in detail the application of GUC to the 
Michaelis-Menten process. We continue our discussion of the No-Pumping Fluctuation Theorem (NPFT) in Section 
[C| first commenting in Section |C 1| on the relation between NPFT and the other, previously known, Fluctuation 
Theorems (FT), followed by outlining the possible application of NPFT to turnstiles in Section |C~2} 



Appendix A: Further Discussions of Geometric Universality of Currents (GUC) 
1. Proof of Geometric Universality of Currents (GUC): Eulerian Approach 



Geometric universality of currents can be also demonstrated using the operator approach. The latter is implemented 
to study the generating function, 



Z(q;T) = 5>(J|T) [] qt , 



(Al) 



depending on the multi- variant argument q = (qi k \{i, k) e Si)- We will focus on the analysis of the T ^> r asymptote, 
where Z(q;T) ~ e TA ( q \ Since §(j), with j = J/T, is related to A(q) by the standard Legendre transform 



S(j)= E 3iklnq* k -A(q*), 
V(i,fc)eg i: j ik =q* k d qik A(q% 



(A2) 



the universality of the LD function, S(j), under the considered transformations is equivalent to the universality of 
A(g). 

Within the respective operator, so-called Doi-Peliti, formalism [3S1I37] (see also [3HHH]) the Master Equation [see 
Eq. (1) of the main body] becomes, dP(n;t)/dt = H(t)P(n;t), where the time-dependent operator H(t) is 



at 



i,k^0 

H(t) = Yl A 0i (af-1)+ 

(0,i)S9i (m)£9i 
(i,0)£9i fceSo 



kit) 



(A3) 



with a\P{n) = P(n_j) and 



b k (t)P(n) = (rij + l)K k (t,n +k )P(n +k ). 



(A4) 



In deriving the last term in Eq. (A3 1 we used ^[. l ' fc ' )eSl di k = 1. Note that H(t) has a very simple form, being 
expressed via the standard "creation" operators a k) which are local and time- independent, and the non-standard 
"annihilation" operators bk(t), which are non-local and depend on time explicitly. The complexity of the system 
stochastic dynamics is hidden in the time-dependent, nonlocal, and nonlinear nature of the "annihilation" operators. 
As we will see, the function, A(q), is not sensitive to the specific form of b k (t), which provides the Euler (operator) 
view of the GUC. 

Extending the standard operator approach developed for generating functions in [34] to the most general case, 
where the rates depend on the occupation numbers and time arbitrarily, we derive 



Z(q-T) = U q P (n), U q {T) = Texp uf dtH q {t)^j 



(A5) 
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Here, Po(n) is the initial distribution, Texp stands for the time-ordered exponential, and the twisted master operator 
H q (t) is 

H q (t) = ]T A o s :(<Zo4-l)+E 5 ^9)&,W 
(CM)6Si i&So 
(i,fc)eSi 

= E 1ikQika\ - at + g io 0io- 



(A6) 



In the long-time limit the spectral decomposition of Eq. (A5) is dominated by the eigenvalue of the time-ordered 



exponential operator with the highest absolute value. Solving the eigenvalue problem U q \^ q ) — z q \^ q ) in the most 
general case does not appear feasible, thus reflecting the complexity of the most general stochastic dynamics of 
interacting particles. However, the non-trivial simplicity in the system that leads to the universal property of the 
generated currents shows itself when we study the eigenvalue problem for the bra-, rather than ket-, eigenstate: 
{^q\Uq = {^q\ z q ■ The bra-eigenvalue problem for the highest absolute value can be treated explicitly by observing 
that for given q the family of operators H q (t), parameterized by time, all share the same bra-eigenstate with the same 



eigenvalue, i.e., (^/ q \H q (t) 



(ty q \i<j q , so that z q 



e w « , which yields A(q) 



The bra-eigenstate is represented by a coherent state (C| = (0| IljeSo e jaj , with (0| denoting the state with no 
particles. The coherent states are parameterized by £ = € So) an d satisfy the following important property: 



(C| a] = (CIO- Due to the latter property a coherent state (C| is an eigenstate of the first term on the rhs of Eq. ( A6 1. 

Therefore, to ensure that (C| is an eigenstate of H q (t) it is sufficient to verify that (C|ct = for all i. Consequently, 
this condition results in a system of linear equations 



ifcfefcCfc — Ci) + 9io{qiO ~ Ci) — 0, 

q = ^oi(qoiCi — 1)7 



MO 



(A7) 



whose solution fully identifies A(q) = uj q 
time-independent parameters, Oik and Aq^ 



Finally, the GUC follows from the fact that (A7) depends only on the 



2. Proof of GUC for Dynamic Pumping 



It appears that the Eulerian derivation of the main GUC presented above allows for a straightforward and very 
simple generalization proving GUC for dynamic pumping. Indeed, the first equation in Eqs. (A7) does not depend 
on the pumping rates Aoi at all. Thus (\P q | = (C|, with C obtained by solving (A7), will also be independent of the 



pumping rates. Therefore, {9 q \ = (CI forms a bra-eigenstate for the time-ordered exponential, defined in (A5) 



(* q \U q (T) = (* q |exp( f dty2\ 0i (t)(q 0i Q - 1). 
Jo 



(A8) 



This results in a universal expression, A(q) = uj q — ^\ XoiiqaiQ — 1), for the general case of time-dependent pumping 
rates. Applying the Legendre transform (A2| to A(q), one arrives at the GUC for dynamic pumping. | 



Appendix B: Counting Statistics in the Michaelis-Menten Process 



As stated in the main body, the counting statistics of the enzymatic reaction events does not provide full quantitative 
information on all kinetic rates of elementary reactions involved in the process. Here we illustrate the general principle, 
using the model shown in Fig. 2 of the main text. We are interes ted in the counting statistics of transitions through 
the link (2, l) 2 only, thus setting all the generating factors in Eq. (A7), except (fci 1 ; to one. To simplify the notations, 
in this paragraph we will be using q instead of <72i 2 - Ovrr focus is on the analysis of A(g) = fc+(Ci — l)j where Ci 



satisfies the following equation, which is a version of ( A7 ) 



*l[S](C2 

fc-i(0- 



- Ci) - k E d 

C 2 ) + fc 2 (gCi 



: -k E> 
C2) - kc£% 



-ke- 



rn 
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This results in the following expression for the leading exponent in the generating function, 

\{q\ = (g - l)k 1 [S\k 2 k+ 

W k 1 [S}{k 2 (l-q)-k c ) + k E (k 2 + k c + k^ 1 )- { ' 

This expression shows the dependence on the relative probabilities, kc/k 2 , k-\/k 2 and kE/ki[S] only, and it does not 
depend on the rates explicitly. Obviously, a similar statement can be made regarding the LD function of the current, 



defined according to the Legendre transform (A2|. To avoid bulky expressions we are not showing §(j) explicitly, 
but rather present the corresponding expressions for the average current and its variance (these are the two major 
characteristics currently available experimentally, while measuring the higher cumulants still represents a significant 
experimental challenge as it requires higher accuracy) 

(j) = qd q A(q)\ q ^ (B3) 
w k 2 [S}(k + /k_) R k 2 + fc_i 

[S] + K M m ' fci 

i((i-0'» 2 ) = (d q qd q A(q) - (qd q A(q)f) \ q=1 (B4) 



2[S] 2 kjk 



+ 



k 2 _([S}+K M M\ 



2 ' 



In deriving Eq. (B3) we assumed that the enzyme and its complex degrade with the same rate k E = kc = fc_, and 
considered the limit of slow degradation, fc_ <C k-i,k 2 , ki[S]. The combination of parameters, denoted in Eqs. ( B3|B4 ) 



by Kmm j is called the Michaelis-Menten constant. Since the average number of enzymes is given by the ratio fc + /fc_ 
the average current (j) in Eq. ( |B2[ ) is given by the famous Michaelis-Menten law for the average enzyme concentration. 
The effect of the enzyme number fluctuations on the variance of the product creation rate is, however, dramatic. It 
has been shown previously [37] that the noise in an MM reaction with a fixed number of enzymes is suppressed, so 
that the so-called Fano factor / = T((j — (j)) 2 ) / (j) is smaller than unity, thus showing sub-Poisson statistics. In our 
case the Fano factor is given by 

/ = 1+ fc_([S]fl M M)' (B5) 

i.e. the underlying statistics is super-Poissonian, with the Fano factor being much larger than one, as k 2 /k- <C 1. 
From the Michaelis-Menten curve for the average current, presented as a function of [S], one can measure only two 
independent quantities, namely k 2 k + /k- and the MM constant Kmm- Measurements of the Fano factor provide 
additional information on the ratio k 2 /k—, thus allowing the value of fc + to be extracted from a combination of the 
two measurements. However, the other rates k\ and k-\ cannot be determined, even if k_ is known independently 
from some complimentary experiment. It is instructive to compare these no-go results with the calculations of the Fano 
factor for the class of enzymatic mechanisms with a fixed number of enzyme molecules, considered in [25 , where it was 
shown that if the number of enzyme molecules is fixed, measurements of the Fano factor are sufficient to distinguish 
among different enzymatic mechanisms, even if on average they produce identical MM curves. Moreover, according 
to [35], it was possible to determine the values of all the kinetic rates for the transitions among the metastable states 
within these complex reactions. Our result shows that for a wider class of reactions, such measurements of the reaction 
event statistics can be insufficient to fix the complete set of kinetic rates. Hence other statistical characteristics that 
are also sensitive to the temporal, rather than only geometric properties of stochastic trajectories (e.g. representing 
statistics of the time intervals between successive events), would be required. 



Appendix C: Further Discussions of No-Pumping Fluctuation Theorem (NPFT) 

1. Relation between NPFT and previously known FT 

NPFT contains fluctuation theorem in its name because of its close relation to other fluctuation theorems formulated 
in the past for systems with DB. Recently, the statistical theory of fluctuations in mesoscopic systems has experienced 
considerable progress due to the discovery of exact results called fluctuation theorems [THU |49l [50] and nonequilibrium 
work relations [5j l51"H53] . however, unlike NPFT, the previously found exact results describe or follow from the 
relations among standard thermodynamic characteristics such as entropy, work, or information. The NPFT represents 
a fluctuation theorem of a new type. It is related to but different from the recently discussed pumping-restriction 
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theorems and related results stated for the average currents [5~4"H5T] in closed networks. The latter pumping-restriction 
theorems, however, cannot be generalized to current fluctuations, neither can they be derived as a special case of the 
solution of Eqs. (17-18) of the main part of this paper. 

Other, previously known, fluctuation theorems can be derived from the general solution of Eq. (17) from the main 
body. If one starts with Eq. (18) of the main body, divides it by q^, multiplies by = e Mfc , and then sums over k, one 
finds that the rhs of the resulting expression coincides with Eq. (17) of the main body. If we additionally assume that 
qi = e Qi and q^ = e^^ hlk ^^ a ^ ak \ where aj, i = 1, . . . N are arbitrary parameters, then the lhs will be identical 
to zero, since the expression under the Yl^j is odd with respect to a permutation of the indices. Considering the 
definition of A(q), the result is equivalent to the following fluctuation theorem for the vector of the total number J 
of particles, transferred through the links 

| e E, °tiJi+Y, i j[(vi-t*3)-(°ti-°'j)]Jij\ — i (CI) 



(CI ) is a Jackson network generalization of the standard fluctuation theorem called the Jarzynski equality [5, 52, 53], 
(e _vV ) = 1, where W is the total work performed by external forces over the system. The Jarzynski equality is 



obtained from Eq. (CI) by setting a>i — jii, and noticing that W = X^Mi"^ 



2. Application of NPFT to Electronic Turnstiles 



The NPFT leads to important restrictions on the strategies for building an efficient electronic turnstile [30] , defined 
as a circuit of quantum dots driven by periodic gate voltages at zero bias between the leads. At fixed gate voltages, 
this system does not sustain an electric current on average. However, periodic driving may induce currents, whose 
magnitude can be controlled with high precision since it is possible to pump a specified number of electrons per 
cycle of the driving protocol [13] • Our solution describes such periodically driven electronic circuits in the regime 
of sufficiently large temperatures, when the number of available states inside quantum dots is large compared to 
the typical number of electrons. In this regime, quantum coherence effects and exclusion interactions due to Pauli 
principle can be disregarded, and electronic transitions through the tunnel junctions, separating the quantum dots 
from the leads and from each other, are thermally activated. Kinetic rates are controlled by applying voltages to the 
leads or gate voltages to the quantum dots, as shown in Fig. 3. Many-body electron interactions inside the same 
quantum dots, e.g. correspondent to the Coulomb electrostatic interaction, are allowed. Our solution predicts that in 
the thermally dominated regime one cannot induce a directed current among the leads or tunnel junctions connecting 
quantum dots, no matter how gate voltages are modulated. Moreover, currents, integrated over the driving protocol 
period, will show counting statistics that are indistinguishable from that observed in the thermodynamic equilibrium. 
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